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Precision Lattice Calculation of D and Ds decay constants 
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We present a determination of the decay constants of the D and Ds mesons from lattice QCD, each 
with a total error of about 2%, approximately a factor of three better than previous calculations. 
We have been able to achieve this through the use of a highly improved discretization of QCD for 
charm quarks, coupled to gauge configurations generated by the MILC collaboration that include 
the full effect of sea u, d, and s quarks. We have results for a range of u/ d masses down to ms/5 
and three values of the lattice spacing, which allow us to perform accurate continuum and chiral 
extrapolations. We fix the charm quark mass to give the experimental value of the rjc mass, and 
then a stringent test of our approach is the fact that we obtain correct (and accurate) values for 
the mass of the D and Ds mesons. We compare fo and /_d, with fx and and using experiment 
determine corresponding CKM elements with good precision. 



I. INTRODUCTION 

Precision calculations in lattice QCD play a crucial 
role in testing our non-perturbative theoretical tools, 
by comparing the results of the calculation with pre- 
cisely measured quantities. In addition accurate cal- 
culations of non-perturbative QCD quantities are very 
important in the extraction of information from anal- 
ysis of experimental data, for example in the deter- 
mination of the Cabibbo-Kobayashi-Maskawa (CKM) 
matrix elements. 

This is most clearly seen in the case of "gold-plated" 
processes, for example the leptonic decay of D^, Dd, tt 
and K mesons. In this process the corresponding me- 
son, with quark content ah (or ah) annihilates weakly 
into a W boson, with a width given, up to calculated 
electromagnetic corrections [HI^, by: 
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Vab is the corresponding element of the CKM matrix, 
and the decay constant fp parametrizes the amplitude 
for W annihilation. By combining a measurement of 
r with an accurate calculation of fp ^ can be used 
to determine Vah- If Vah is known from elsewhere we 
can use (HJ to get a value for fp. 

The decay constant fp is conventionally defined to 
be a property of the pseudoscalar meson, calculable 
in QCD without QED effects, and is given by: 



(2) 



The calculation of fp is a hard non-perturbative 
problem, which at present can only be done fully with 
lattice QCD. There are very precise experimental mea- 
surements for the leptonic decay rates in the case of 
the TT and K, and new results are appearing for D 
and Ds, which make the calculations a highly non- 
trivial test of lattice QCD, and ultimately of QCD 
itself. This tests are important to give us confidence 
in similar lattice QCD predictions of matrix elements 
in B systems, for which experimental results are much 
harder to obtain. 



II. IMPROVED STAGGERED QUARKS 

We use HISQ staggered quarks in the valence sector, 
whereas the sea quarks are ASQTAD staggered quarks 
with the fourth root trick [1, U [H| . 

Staggered quark actions suffer the doubling prob- 
lem: there are four "tastes" (non-physical flavours) 
of fermions in the spectrum, which couple through 
taste-changing interactions. These are lattice artifacts 
of order a^, involving at leading order the exchange 
of a gluon of momentum q « tt/o. Although quite 
large in the original one-link (Kogut-Susskind) stag- 
gered action, such interactions are perturbative for 
typical values of the lattice spacing, and can be cor- 
rected systematically a la Symanzik. By judiciously 
smearing the gauge field we can remove the coupling 
between quarks and high momentum gluons. 

The most widely used improved staggered action 
is called ASQTAD, and removes all tree-level a^ dis- 
cretization errors in the action [1, 0, [1] ■ 

The HISQ (highly improved staggered quarks) stag- 
gered Dirac operator involves two levels of smearing 
with an intermediate projection onto SU{2>). It is de- 
signed so that, as well as eliminating all tree-level a^ 
discretization errors, it further reduces the one-loop 
taste-changing errors (see ^] for a more detailed dis- 
cussion.) This action has been shown to substantially 
reduce the errors associated with the taste-changing 
interactions fol. [lOl. . 

When we put massive quarks on the lattice, the dis- 
cretization errors grow with the quark mass as powers 
of am. Therefore to obtain small errors we would 
need am <C 1. For heavy quarks this would require 
very small lattice spacings. On the other hand, to 
keep our lattice big enough to accommodate the light 
degrees of freedom, we need La ^ m~'^ . The fact 
that we have two very different scales in the problem 
makes difficult a direct solution. What we can do in- 
stead is to take advantage of the fact that m is large, 
by using an effective field theory (NRQCD, HQET). 
This program has been very successful for b quarks 
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The charm quark is in between the hght and heavy 
mass regime. It is quite hght for an easy apphcation 
of NRQCD, but quite large for the usual relativistic 

quark actions, arric ^ 1. However, if we use a very 
accurate action (HISQ) and fine enough lattices (fine 
MILC ensembles) , it is possible to get results accurate 
at the few percent level. A non-relativistic analysis 
shows that for HISQ charm quarks the largest remain- 
ing source of error is due to the quark's energy, and 
can be further suppressed by powers of v/c, where v 
is the typical velocity of the quark in the system of 
interest, simply by retuning the overall coefficient of 
a term called Naik term to impose the correct rela- 
tivistic dispersion relation c^(p) = 1 for low lattice 
momemtum p. 

One advantage of the use of a relativistic action is 
the existence of a partially conserved current, which 
implies the non-renormalization of the lattice result 
for fp. We can extract fp from the PC AC relation 
for zero momentum meson P: 

fpml = {ma + mt) (0|a75&|P) (3) 





III. RESULTS 




Lattice/sea 


valence 


ri/a 


uoanii, uoanis 


ami, ams, anic, 1 + e 




16^ X 48 






0.0194, 0.0484 


0.0264, 0.066, 0.85, 0.66 


2.129(11) 


0.0097, 0.0484 


0.0132, 0.066, 0.85, 0.66 


2.133(11) 


20^ X 64 






0.02, 0.05 


0.0278, 0.0525, 0.648, 0.79 


2.650(8) 


0.01, 0.05 


0.01365, 0.0546, 0.66, 0.79 


2.610(12) 


24^ X 64 






0.005, 0.05 


0.0067, 0.0537, 0.65, 0.79 


2.632(13) 


28^ X 96 






0.0124, 0.031 


0.01635, 0.03635, 0.427, 0.885 


3.711(13) 


0.0062, 0.031 


0.00705, 0.0366, 0.43, 0.885 


3.684(12) 



TABLE I: MILC configurations and mass parameters used 
for this analysis. The 16'^ x 48 lattices are 'very coarse', 
the 20^ X 64 and the 24^ x 64, 'coarse' and the 28^ x 96, 
'fine'. The sea asqtad quark masses (/ = u/d) are given 
in the MILC convention where uo is the plaquette tad- 
pole parameter. Note that the sea s quark masses on 
fine and coarse lattices are above the subsequently deter- 
mined physical value [TJ- The lattice spacing values in 
units of ri after 'smoothing' are given in the rightmost 
column [l6l. Hsj. The third column gives the HISQ valence 
u/d, s and ^masses along with the coefficient of the Naik 
term, 1 -I- e, used for ^quarks 

We use gluon field configurations including 2 + 1 
flavours of sea quarks generated by the MILC collab- 
oration [l^, [l^ . The parameters of the ensembles 
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FIG. 1: Charmonium spectrum obtained with HISQ 
charm quarks on fine MILC lattices (blue crosses with er- 
ror bars) against experiment (black lines). 



we have used for both the sea and the valence sec- 
tors are in table [D The lattice results are converted 
to physical units through the heavy quark potential 
parameter ri, as determined by the MILC collabora- 
tion (table m fl^). The physical value of ri is de- 
termined from the T spectrum calculated in NRQCD 
with b quarks on the same MILC ensembles [13], with 
the result ri = 0.321(5) fm, rf^ = 0.615(10) GeV. 

We use multiple precessing random wall sources, 
which gives a 3-4-fold reduction in statistical errors 
with respect to conventional local sources. 

The mass of the charm quark is fixed by adjust- 
ing the mass of the "goldstone " r]c to its experimen- 
tal value. The hght {u/d) and strange quark masses 
are fixed using the experimental values for the masses 
of TT and K. Our results use masses for the u and 
d quarks that are substantially larger (by a factor of 
around three) than the real ones. In order to get phys- 
ical answers we extrapolate to the correct u/d mass 
using chiral perturbation theory. Once the masses 
have been thus fixed, there is no remaining freedom to 
change any parameters, and in particular the results 
we obtain for the masses of heavy-light mesons are 
a stringent test of our method. In figure [T] we show 
the spectrum of charmonium. We obtain an hyperfine 
splitting of 111(5) MeV (experiment 117(1) MeV) We 
have made no attemp as yet to optimize the calcula- 
tion of the excited states. 

In addition to the chiral extrapolation, we have sys - 
tematic errors coming from a variety of sources [19|, 
among them from the finite lattice spacing. Because 
we have three different lattice spacings and very pre- 
cise data, we can extrapolate to the continuum limit. 
This extrapolation is linked to the chiral extrapola- 
tion through discretization errors in the light quark ac- 
tion. We therefore perform a simultaneous bayesian fit 
for both chiral and continuum extrapolations, allow- 
ing for expected functional forms in both. We tested 
the validity of the method by fitting hundreds of fake 
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FIG. 2: Masses of the and Ds meson as a function of 
the u/d quark mass in units of the s quark mass at three 
different values of the lattice spacing. The very coarse 
results are above the coarse and the fine are the lowest. 
The lines give the simultaneous chiral fits and the dashed 
line the continuum extrapolation as described in the text. 
Our final error bars, including the overall scale uncertainty, 
are given by the shaded bands. These are offset from the 
dashed lines by an estimate of electromagnetic, niu 7^ rnd 
and other systematic corrections to the masses. The ex- 
perimental results are marked at the physical md/rns- 



datasets generated using staggered chiral perturbation 
theory with random coupUngs. We fit simultaneously 
to the masses and the decay constants, that is, we fit 
TO^, rriK, /it and /k simultaneously, and similarly for 
mo, toDs , Id and /^^ . We present some of the results 
in figures [2] and [31 

We get an excellent agreement with experiment 
for the masses: m/j^ = 1.963(5) GeV (experiment 
1.968 GeV), and tod = 1.869(6) GeV (experiment 
1.869 GeV). Our calculation also reproduces correctly 
the difference in binding energies between a heavy- 
heavy (rjc) and a heavy-light {mo and mo^) state: 
{2mr)^ — rriri^) / {2mr) — ti^^) = 1.249(14) (experiment 
1.260(2)). Our charm quark action is the first one to 
be accurate enough to do this calculation (which also 
cannot be done, for example, in potential models.) 

We also have agreement with experiment for the 
light-light decay constants [l§|. The result for the 
ratio is very accurate, Jk/ f-n = 1.189(7), and shows 
tiny discretization effects (figure [3]) . Combining this 
ratio with experimental leptonic branching fractions 
[UllOl we get Vus = 0.2262(13)(4), where the first 
error is theoretical and the second experimental. This 
gives the unitarity relation 1 — V^^ — V^^ — V^j, = 
0.0006(8). 

Our results for the heavy-light decay constants 
are 4-5 times more accurate than previous lattice 
QCD results and existing experimental measurements: 
fo, = 241(3) MeV, fo = 208(4) MeV, and a ratio of 
foj fo = 1.162(9) (see figure[5]). For the double ratio 
ifoj fD)/{fK/ fyr), which is estimated to be close to 
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FIG. 3: Results for the D, Ds, K and tt decay constants 
on very coarse, coarse and fine ensembles, as a function of 
the u/d quark mass in units of the s quark mass. The chiral 
fits are performed simultaneously with those of the corre- 
sponding meson masses, and the resulting continuum ex- 
trapolation curve is given by the dashed line. The shaded 
band gives our final result. At the left are experimental 
results from CLEO-c [H, ^23] (on the left with the r decay 
result above the fi decay result for Ds) and BaBar [S^l (-Ds 
only) and from the Particle Data Talales for K and vr. 
For the K we have updated the result quoted by the PDG 
to be consistent with their quoted value of Vus ■ 
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FIG. 4: Ratio of decay constants fK/fn on very coarse, 
coarse and fine ensembles, as a function of the u, d quark 
mass in units of the s quark mass. 
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FIG. 5: Ratio of heavy-light decay constants fosl fo on 
very coarse, coarse and fine ensembles, as a function of the 
u, d quark mass in units of the s quark mass. 



1 from low order chiral perturbation theory [2l|, we 
get a value of 0.977(10). 

The experimental leptonic branching rates, together 
with CKM matrix elements determined from other 
processes (assuming Vcs = Vud) give fos — 264(17) 
MeV for udecays and 310(26) MeV for r decay from 
CLEO-c [22^ and 283(23) MeV from BaBar (H, and 
for fo 223(17) MeV from CLEO-c for ^ decay [H]. 
Using our results for fo^ and fp/fp and the exper- 
imental values from CLEO-c [22] for /i decay (since 
the electromagnetic corrections are well-known in that 
case) we can directly determine the corresponding 
CKM elements: V^s = 1-07(1)(7) and Ks/Kd = 
4.42(4)(41). The first error is theoretical and the sec- 
ond experimental. The result for Vcs improves on the 
direct determination of 0.96(9) given in the Particle 
Data Tables [1]. 

Our calculation is precise enough that we can see 
the difference between mB^(m;) — mB{mi) in the bot- 
tom sector and the similar quantity in the charm sec- 
tor mp^{mi) — mp{mi) (figure [6]). These mass differ- 
ences are small compared to the absolute masses of 
the states, and should be the same in the infinitely 
heavy quark limit. We can see that our calculation 
correctly reproduces the small difference due to the 
finite value of the mass of the charm and bottom. 
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FIG. 6: Binding energy difi^erences of heavy-light mesons 
with strange valence quarks and those with u, d va- 
lence quarks, for states with a ^quark and systems with 
a h quark, — rno and — rns- Data are shown 
for the coarse and fine ensembles in the case of D and Ds, 
and for the coarse for the B and Bs. 



controlled joint chiral and continuum extrapolation. 

We can calculate accurately the mass of heavy-light 
systems, which provide a stringent test of the calcu- 
lation. We can calculate precise values for the decay 
constants of pseudoscalar heavy-light mesons (as well 
as light-light mesons), and especially for the ratio of 
such decay constants. 

The very precise calculation of the masses of heavy- 
heavy pseudoscalar mesons should make possible a di- 
rect lattice determination of the mass of the charm 
quark. Because we use the same relativistic action 
through the calculation for both the charm and the 
light quarks, we can also obtain a very precise value 
for the ratio mc/rris, and therefore if rric is deter- 
mined through another method use the ratio to get 
rris- We are also working (in collaboration with the 
Karlsruhe group) on a new method for the determi- 
nation of rric by combining continuum perturbation 
results with lattice data. 

Another quantity which we plan to calculate in the 
near future is the leptonic decay width Tf.+e-i''P)j as 
well as the semileptonic form factors for D — > nlv, 
D Klv. 



IV. CONCLUSIONS AND OUTLOOK 

We have shown that the use of a highly improved 
relativistic action on fine enough lattices is capable of 
delivering very precise results on systems with a charm 
quark. The high statistical accuracy of our data com- 
bined with calculations at several values of the lattice 
spacing and light quark masses allows us to make a 
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